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ABS'^5ACT 

This module is part of a series designed to be used 
bv life science students for instruction in the application of 
physical theory to ecosystem operation. Most modules contain computer 
proarams which are built around a particular application of a 
physical process^ This module describes heat transfer processes 
involved in the exchange of heat between an organism and its 
environment. Emphasis is placed on conduction, convection, and 
evaporation, three of the physical processes that affect the 
physiological, behavioral, and ecological activities of all 
organisms. Padiation transfer is described in a separate aodule. Four 
examples taken from the ecoloaicai literature are used to emphasize 
the bioloaical effects of temperature and energy budgets. Three 
systems are used to illustrate heat transfer principles and their 
biological importance. These include heat flow in soil, a leaf, and a 
lizard. Ecological examples are used throughout the text and problem 
set. Calculus and first-year physics are prereguisites for the 
module. (Aut hor/CS) 
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PREFACE 

The way in which an organism exchanges heat with its environment 
explains several characteristics of its behavior and its preferred habitat 
This module presents a thorough introduction to heat transfer processes 
and assumes the reader has a background in calculus and first-year physics 
Emphasis is placed on conduction, convection and evaporation. Radiation 
transfer is described in its own module. Ecological examples are used 
throughout the text and the problem set. 
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INTRODUCTION 

The purpose of this module is to further elucidate the physics of 
; the heat transfer processes: radiation, evaporation, conduction, and 
r convection. An understanding of these processes and their Interactions 
will provide a clearer ecological interpretation of the thermal energy 
environment. Two ideas need to be recalled about the First Law of 
Thermodynamics (Zemansky and Van Ness 1966, Stevenson 1979a): first, that 
the heat energy budget is based on the conservation principle of the First 
Law, and second, that this law can be applied to any system of arbitrary 
boundaries. ; Although organisms seem a natural choice, Collins et al. (1971) 
have found tfeis^ nasal passageway to be a useful system for study of nasal 
vapor recovery. In a previous module (Stevenson 1979Z?) , we considered 
a stream, a: leaf, and a spider as well-defined thermodynamic systems for 
understanding the biological importance of heat fluxes. 

Before we begin to describe the mechanisms of heat transfer, four 
examples are taken from the ecological literature to emphasize the 
biological effects of temperature and energy budgets. As the examples are 
presented, the reader should ask, "How does the thermal environment 
influence the distribution, abundance, and life history strategies of the 
organism being described?" 

Boylen and Brock (1973) have studied the benthic algae of the 
Firehole River in Yellowstone National Park. They were interested in 
documenting the effect? of the heated water entering the river from 
geysers. The Firehole was coldest during June when the discharge was 
largest from the snowmelt. Diatoms in the genera Nitzchia^ Snydedia^ 
Rhopalodiaj Cocconeis^ and Comphonema were present at the sampling 
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stations with cooler water temperatures, while green algae, Spirogyra^ 
Oedogoniurrij Cladophora^ and S tig eoa Ionium ^ were dominant at the warmer 
stations. Biomass was 15 times larger and growth rates five ttmes greater 
tin the heated than the unheated section of the river, Boylen and Brock 
^eel that the Firehole River represents a unique opportunity 
to study the biological consequences of elevated water temperature which 
are similar to the increases caused by power plants. 

The acquisition of nutrients is of fundamental importance to plants. 
Chapman (197A) studied the absorption of phosphate along thermal gradients, 
where: 1) the average soil temperature changed and the variability was 
constant; and 2) where the mean was the same hut the fluctuation or 
variability was changing. The maximum rate of phosphate absorption was 
found to correlate well with the mean soil temperature for each species. 
He found that cold-adapted species also increased the root-to-shoot ratio 
presumably in order to compensate for the slower absorption rates at colder 
temperatures. Species from fluctuating environments showed a greater rate 
of acclimation to phosphate absorption than did species from more constant 
environments. These results suggest that plants from different 
environments try to maintain similar phosphate uptake rates, which is 
accomplished with different physiological adaptations, and that plants can 
partition their enerc;y resources (root and shoot biomass) to achieve this 
balance . 

Our third example concerns temperature adaptations in amphibians. 
Snyder and Weathers (1975) hypothesized that the variability tolerated in 
body temperatures of amphibian species would be correlated with the 
variability in air temperaturjes . To test this idea, the difference 
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between the Tnaximum and lower lethal body temperature for 11 species was 
taken from the literature and plotted against the difference in the high 
-and low mean monthly air temperatures for a 10-yr period. Figure 1 
presents their results and is consistent with their hypothesis. 

Several authors (Schneirla et al. 1954, Scherba 1960, Jackson 1957) 
have investigated the effects of microclimate on ant activity and nest 
structure. Nest and bivouac temperatures are important for egg, larva 
and pupa development, Golley and Gentry (1964) found that soil surface 
. temperature was an important factor influencing the activity of the 
Florida harvester ants. Figure 2 shows that Pogonomyrmex badius was not 
active in June and July when the soil surface became very hot during 
midday hours. Also there was little activity in February, when soil 
surface temperatures were cold. Carlson and Gentry (1973) further found 
that continual shading of this species' entry holes causes a large 
increase in the migration of colonies to new nest sites. Sunshine is 
probably an important energy contribution needed to raise the ant's body 
and gallery temperatures. 

It is sometimes useful to take a more general outlook. In that 
spirit, here are six ecological questions which could be more easily 
interpreted and answered with a sound knowledge of heat transfer physics: 

1) Does the climate limit the geographical distribution of an 
organism? 

2) What selective advantage is there to being a particular color, 
shape, or size? 

3) What microhabitats are available throughout the day "to maintain 
a thermal balance? 
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ENVIRONMENTAL TEMPERATURE VARIATION 

Figure 1 , The relation of increase in range of temperature tolerance and 
increase in environmental temperature variation (see text for methods of 
determining these values). The continuous line is a line of equality relat- 
ing increase in temperature tolerance and environmental temperature variation. 
The dashed line has been fitted to the data with the method of least squares. 
The numbers indicate the species included. These are: 1) Bufo alvarius^ 2) 
B. debilis^ 3) B. marinus^ 4)S. marmorcus^ 5) B. mazatlanensis ^ 6) Hyla 
Tegillay 1)H. smithi^ 8) Leptodaatylus melanolus^ 9)Ptemohyla fodiens^ 
10) Scaphiopus haxmondi^ 11) Smilisca baudini. (From Snyder, G.K. and W.W. 
Weathers, 1975, p. 98). 
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■^^g"^^ ^- Surface activity of Pogonomyrrnex badius in relation to temperature 
of^'jthe ground surface and time of day. Active ants on the mound '.are indicated 
by;crosses; inactive mounds are indicated by open circles; half-closed circles 
denote very slight activity. The months: represented , beginning with the high- 
es,l:..|.curve at 12 noon, are June, July, May, October and February. (From Gollev 
F.Bi and J.B. Gentry, 1964, p. 224.) 



4) Is there a daily activity pattern that limits prey-predator 
interactioris or foraging time? 

> 5) What selective advantage is there to be a poikilotherm or 
homeotherm? 

6) How does the thermal balance affect photosynthetic and 
transpiration rate of a plant? 

The reader should be cautioned that there are other ecological 
factors such as competition^^ predation, and coevolution, which may be 
greater selective influences and which may interact with selective 
pressures of the physical environment. Predictions from models incorpora 
ting physical processes can often be tested quantitatively which allows 
the investigator to compare the importance of the physical process under 
consideration with other ecological factors. 

nzAT train;sfer processes 

Radiation 

Radiation is a process of energy transfer that requires no interven- 
ing medium. The dual nature of radiation, its particle (photon) and 
wavjslike properties have important biological consequences. All matter . 
radiates energy as individual photons (Kreith 1973) but for heat transfer 
problems the wavelike description is more useful. This is because 
energy at different wavelengths internees with matter in different ways 
(but see Problem 2). Figure 3 fihows the different wavelengths of electro 
magnetic spectrum. 

Since radiation travels at the speed of light c^, the product of 
the frequency f and thie wavelength A is a constant: 
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Figure 3, The electromagnetic spectrum. (From Byers, H.R., 1974, p, 




Figure 4. Schematic representation of blackbody emission spectra, 
as a function of wavelength, for temperatures of 6000 
and 270 K. (From Lowry, W.P., 1969, p. 17.) 
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c, = fX. (1) 

All objects, living and nonliving, radiate thermal energy. The amount and 
.kind (wavelength in the electromagnetic spectrum) of energy depend on the 
temperature and physical characteristics of the radiating body. Figure 4 
shows the energy radiated by the sun and the earth as a function of 
wavelength* It should be clear from the figure that the portion of the 
electromagnetic spectrum radiating from the sun occurs at much shorter 
wavelengths than that of the earth. In fact, the sun's peak radiation on a 
wavelength plot is in the green part of the visible spectrum, while the 
earth's radiation is completely in the infrared region. 

Three physical laws are associated with Fig. A. Planck's Law gives 
the energy emitted, , as a function of the wavelength, X, if the 
temperature of the radiating body is known: 

= ciX-5 [exp(c2/XT) - l]"-^ (2) 

where 

is the amount of energy emitted in the band A to A + dX (J m*"-^),^ 
T is the blackbody temperature (K), 
X is thf* wavelength (m), 
c^ is 27ThCL^, 
C2 is hcL/K-b, 
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where 

h is Planck's constant = 6.63 x 10"34 (j g), 

is Boltzmann's constant = 1,38 x 10"'23 (j k"^), 

and 
#^ 

Cl is the speed of light = 3*00 x 10^ (m s""^). 

It is obvious that no one has actually measured the surface temperature of 
the sun. In this instance, the sun's temperature was deduced from its 
electromagnetic spectrum. Second, from Planck's Law, it can be 
demonstrated (Robinson 1966 or Roseman 1978) that the wavelength at the 
maximum radiated energy is only a function of the temperature of the 
radiating body, 

\ 2>897 ,^-3 

A - X 10 , 

max T (3) 

where 

"^max ~ wavelength of maximum radiation (m), 

and 

T = temperature (K), 

Third, the S tef an-Boltzmann Law provides the even more remarkable result 
that the total energy emitted, Qg , by a radiating body (the area under the 
curve (Fig. 4)) is proportional to the fourth power of the surface 
temperature Tg . 

Qe = £^(273 + Tc)^, (A) 
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where 



is the radiant energy emitted from the surface (W ni"2)^ 
e is the emissivity (range zero to one), 

a is the Stephan-Boltzmann constant (5,67 x 10*^8 y j^-A^ 
8.17 X 10-^1 cal min"^cni"2K-^), 

and 

Tg is the surface temperature (°C). 

Roseman (1978) gives a more complete discussion of these laws. 

Kirchoff studied the emissivity and absoiptivit;- of materials. He 
showed that the energy absorbed at any specific wavelength X at constant 
temperature is equal to the energy emitted at X. That is, if the 
absorptivity of a surface a(X) represents the fraction of incident 
radiation absorped at X and e(X) is the emissivity, the radiation emitted at ^ 
divided by the radiation emitted from a blackbody, then a(X) = e(X) for each 
wavelength. If e = 1 for all wavelengths, then the object is said to be a 
blackbody. Although no system is a perfect blackbody, many are 
approximately so. In the infrared region, 3 to 100 ym, most objects such 
as vegetation, soil and water behave like blackbodies (e = !)• Table 1 
gives the emissivities of some plants and animals. 

This, however, is not true in the visible part of the spectruir. In 
fe-iral, any wavelength can be absorbed, reflected, or transmitted. This 
is written: 

f- 

a(X) + r(X) + t(X)>''l •• ■ (5) 

t ■ 
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Table 1. Long Wave Emissivities (from Monte ith, J.L., 1973, p. 68.) 







LEAVES 






Species 




Average 





Maize (Zca rr,c'js) 
Tobacco CUcoiicKa tabacur:) 
Snap bean (Pksseclus vulgaris) 
Cotton iCcssypiun hirsutur. Deltapine) 
Sugar cane (Sac^hzjrur, officinar-un) 
Poplar (Populus fre-nontii) 
Geranium ('Fclci:rgoniun dcmeszicurJ 
Cactus (OpUKtia rufidz) 



94.4 
97.2 
93.8 
96.4 
99.5 
97.7 
99.2 
97.7 



0.4 
0.6 
0.8 
0.7 
0.4 
0.4 
0.2 
0.2 



ANIMALS 



Species Dorsal Ventral Average 



Red squirrel (Tanicispiurus 95-98 97-100 

h'udsor.icus ) 

Gray squirrel (Soinrus aaroZiKeKsis) 99 99 

Mole (Sccilopu.- aquaziauc>) 97 

Oeer mouse (^erz^rjscus sp.) — 94 

Grey wolf 99 

Caribou TOO 

Snov/shoe hare 99 

Man (Horr^o sapieyis) 98 



Table 2. Reflectivities of Biological Materials for Solar Radiation 
(from Monteith, J.L., 1973, pp. 66-67.) 



LEAVES 

Reflection coefficients r [%) for Solar Radiation 

Species Upper Lower Average 



Mai.-e (Ze.2 -r-:.€'' 29 

Tobacco {:n-j:!tiaKa tc::-^:irj 29 

Cucumber CCu::uric s::t:v::, 31 

Tomato Clucczycrcz-jsK ecj^uZeKtu-^J 28 

Birch (Setu^d aZb::) 30 33 32 

Aspen fP::puius trc-^Zc'dcs: 32 36 34 

Oak (Qucrcuc alb::; 28 3-1 30 

Elm (^Jlm^s '"ra) .-j 24 31 28* 

. ^ y, 
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Table 2. (Cont'd.) 



VEGETJ>TION — MAXIMUM GROUND COVER 



Farm Crops 


Daily Mean 


Farm Crops 


Daily Mean 


Grass 


24 


Wheat 


22 


Sugar beet 


26 


Pasture 


25 


Barley 


23 


Ba rl ey 


26 


Wheat 


26 


Pineapple 


15 


Beans 


24 


Sorghum 


20 


Ma i ze 


18 - 22 


Sugar cane 


15 


Tobacco 


19 - 24 


Cotton 


21 


Cucumber 


26 


Groundnuts 


17 


Tomato 


23 






Natural Vegetation 


Daily Mean 


Natural Vegetation 


Daily Mean 


Heather 


14 


Natural pasture 


25 


Bracken 


24 


Derived savanna 


15 


Gorse 


18 


Guinea savanna 


19 


Maquis, evergreen scrub 


21 






Forests and Orchards 


Daily Mean 


Forests and Orchards 


Daily Mean 


Deciduous woodland 


18 


Eucalyptus 


19 


Coniferous woodland 


16 


Tropical rainforest 


13 


Orange orchard 


16 


Swamp forest 


12 


Aleppo pine 


17 





ANIMAL COATS 



Mamma 1 s 


Dorsal 


Ventral 


Average 


Red squirrel (Tcniasciurus hudsonicus) 


27 


22 


25 


Gray Squirrel (Sciui^-us carolincncds) 


22 


39 


31 


Field mouse (riicrotus per.KsitlvanicusJ 


11 


17 


14 


Shrew (Sorex sp, ) 


19 


26 


23 


Mole (Scalenus ciquc.zicus) 


19 


19 


19 


Graj' fox (Vrocjor, cine^'^i^o crger.tus) 






34 


Zulu cattle 






51 


Red Sussex cattle 






17 


Aberdeen Angus cattle 






11 


Sheer) ••'oatherod '''leocG 






26 


Newly shorn fleece 






42 


Man (Mono sayieKc) Eurasian 






35 


Negroid 






IS 


Birds 


Wing 


Brtast 


Average 


Ca rd i na 1 ( Eich-'.oKCKKLi ccrdi riUs .■ 


23 


40 




Bluebird 


27 


34 




Tree swallow 


24 


57 




Magpie 


19 


46 




Canada goose 


15 


35 




Mallard duck 


24 


36 




Mourning dove 


30 


39 




Starling (Stun-^r, t^c/^.u'-i"-; 






34 


Gloucous-winged gull (:.r."i<3 jl-ur-.^c^'t-nr, ) 






52 
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where 

a(X) = absorptivity at wavelength X, 
r(A) = reflectivity at wavelength A, 
t(X) = transitivity at wavelength A, 

and each term is between zero and one (see Siegel and Howell 1972 for a 

more complete discussion). Table 2 gives some values for r averaged over 

the entire spectrum for different ecological systems. Values of a, r, t 

are also given by Porter (1967) for animals and Monteith (1973) for 

plants and animals. 

As an example of radiation flux, let us calculate the radiant energy 

emitted by the cactus, Opuntia rufida.. To do this, we must specify the 

emissivity £ and the surface temperature Tg. In Table 1, £ = .977 for 

0. rufida and, if we take T =10 ^C, then: 

s 

Qe = .977 X 5.67 x 10-8 (iq + 273)^ = 355.5 W m"2. (6) 

If we wish to know the total heat loss per second and (1 watt = 1 J s ^) 
the plant, we must multiply our answer by the surface area of the cactus. 

The calculation of absorbed radiation 0^, although straightforward 
with some assumptions, is a lengthy task because several longwave and 
shortwave components must be included. Gates and Stevenson (1979) have 
devoted an entire module to this topic to which the interested reader 
should refer for additional information. 
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Conduction 

Conduction, G, describes the physical process of molecular thermal 
tenergy flow within a solid, fluid or gas. In fluids and especially gases, 
potion of the medium usually makes the process of convection a more 
important mechanism of heat transfer than conduction. Heat flows by 
conduction when nearby molecules have more internal energy (higher temperature) 
and thus a greater mean kinetic energy. Energy can be transferred by 
molecular collisions (fluids) or by diffusion (solids) (Kreith 1973, pages 
4-5). 

In conduction, the rate of heat flow depends on the thermal 
conductivity of the material, the area through which the heat flows and 
the temperature gradient in the material. Equation (7) describes the 
relation as : 

G = -kAef (7) 

where 

G is conduction (W), 

k is thermal conductivity (W ra"-^ °C"M, 

is the area through which the heat is flowing (m^), 

and 

dT 

dx is the thermal gradient (°C m"^). 

The negative sign in equation (7) is to indicate that the direction of 
heat flow is from regions of higher temperature to lower tempejrature 
(opposite the temperature gradient). This is shown in Fig. 5. 
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Oir*<i*OA ol Hmi Flow OwtctKM •( H»<if Flow 








Fig. 5. Sketch illustrating sign convention 

for conduction heat flow. From Kreith, F. 
1973. P. 8. 

The thermal conductivity depends on the molecular structure of the 
material. Some cooking pots have copper-coated bottoms because this 
metal transfers heat better than other materials. Home insulation 
materials reduce the flow of heat because they trap air which has a low 
thermal conductivity. The hollow centers of the hair shafts of some 
mammals reduces the conduction of heat along the fiber. The thermal 
conductivity can be calculated by measuring the heat flow when the 
temperature slope (Fig. 5) is one. 

As an example of conduction, we will examine the heat flow between 

the alligator, Alligator mississippiensisj and substratum. We assume that 

the skin is Lhe same temperature as the ground and that the animal is in 

good thermal contact with the surface (there are no air pockets reducing 

heat flow). The conduction exchange can now be calculated as follows: 

The thermal conductivity of fat is 0.20 W m"^ °C"^ while soil-rock range 

of conductivity is between 0.556 and 3.25 W m"^ ^C"^. Since the fat layer 

of 0.7 cm has a lower conductivity, it controls the rate of heat flow. 

2 

Here we also assume that 0.5^ of the alligator at T^ - 22 °C are in 
contact with the ground at T^ = 30 °C. In steady state, the conductive 
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heat flux is then a linear function of the potential, as given in equation 

(7). . 

kA 

• ° = -ir^h - \^ = •'^O.OO?'' (22-30) = -114 W, - (8) 

-where 

G is conduction (W, negative since energy is flowing from the ground 
to the animal), 

is body temperature of the alligator (°C), 
k is the thermal conductivity of fat (W m"*^ °C~"^), 
T„ is ground temperature (^C), 

and 

Ax is thickness of layer (m) • 

Before proceeding to a discussion of convection, it is important to 
talk about an alternative view of conduction. It is possible to rewrite 
equation (8) as 

G = K (T^ - T ). 

g b g (9) 

The rate of heat flow is then equal to a conductance times a potential, 
the temperature difference. The conductivity, area of contact and length 
have all been subsumed into Kg. When working heat transfer problems, it 
is also common to speak and think of resistances to heat transfer. A 
resistance is simply the reciprocal of the conductance. If there is a 
large resistance to heat flow, there is a low conductance. Monteith 
(1973) and Campbell (1977) have adopted the resistance approach for 
describing mass and momentum fluxes as well as heat transport.* This 
concept originated with Ohm's Law where the electrical current flow is 



30 



16 



equal to the potential or voltage drop divided by the resistance of the 
material t Thus, the correspondence of a flux driven by a potential is 
called Ohm's Law analogy. Kreith (1973, chapters 1 and 2) shoxk how 
eradiation, convection and evaporation can also be thought of in these 
•terras. It should be remembered though that in general the Ohm's Law 
analogy, due to the complications introduced by turbulent mass and 
turbulent energy transfer, is an approximation (see Tennekes and 
Lumley 1973). 

Conduction is often ignored in the heat transfer of plants and 
animals because one must know or raake. assumptions about the temperatures, 
boundary layer thickness, thermal conductivity, and contact resistance. 
There are not many actual measurements reported in the ecological 
literature but Mount (1968), Derby and Gates (1966), Gatesby (1977), and 
Thockelson and Maxwell (1974) offer some interesting examples of the 
importance of conduction. 

Convection 

'^Convection" is the transfer of heat between solids and fluids (i.e., 
gases and liquids) or when fluids of different temperatures are in 
contact. Conduction takes place with nearby particles on the molecular 
level but the additional factor of the circulation of the fluid 
distinguishes convection from conduction. When there is no motion in the 
fluid except that caused by density gradients and by the resultant buoyant 
forces, this process is called *'free convection." If the fluid is moving 
relative to the other fluid or solid, the process is called "forced 
convection." If the flow is turbulent, besides the molecular processes. 
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there will be increased heat transfer caused by the bulk transport of 
fluid. Also, at certain fluid velocities, free and forced convection may 
-►both contribute to the convection term. Indeed convection is comprised of 
several complex physical processes occurring simultaneously • 

Water and air are the two fluids that commonly interest the 
thermobiologist . Most aquatic organisms are at ambient water temperature; 
hence, convective transfer is unimportant because there is no temperature 
difference. Marine mammals, birds, a few large fish and turtles, however, 
maintain body temperatures above water temperature, which requires thick 
insulation (Bartholomew 1977, Schmidt-Nielsen 1975). This is because the 
high specific heat (A. 18 x iO^ J kg"^ OC"^) and thermal conductivity (59.8 
W m"-^ °C"^) of water create a large heat loss. The extra heat transferred 
due to the movement of the fluid has not been examined for many organisms 
(but see Eskine and Spotila 1977 and Lueke et al. 1976). Although the 
specific heat and thermal conductivity of air are much less than for 
water, convection is still a significant mode of heat transfer for 
terrestrial organisms. 

As with conduction, convective heat flow can be thought of as a 
potential (the difference in temperature between the surface of the object 
and the surrounding fluid) times a conductance. The conductance is most 
commonly written as the product of two terms as in equation (10): 

C = he Ac (Ts - T), (10) 

where 

■ V 

C is the convective heat flux (W), 

h^ is the convective heat transfer coef f icient • (W m"^ ^C"l), 
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Ac is the area at object in contact with the fluid (m^), 
Tg is the surface temperature of the object (°C), 

and ^ 
T is the fluid temperature (°C)» 

Usually, A^, Tg and T can be evaluated. (One particular problem, however, 
is the surface temperature of a furry object. Kowalski (1978) has shown 
that in this case it is best to define the system boundary at a solid 
surface where a temperature can be measured.) The heat transfer 
coefficient, on the other hand, is often a difficult parameter to estimate 
in the natural environment because of turbulence (Kowalski and 
Mitchell 1976, and Noble 1975). Generally, it is established as follows: 

he = ^± (11) 
D 

where 

h^ is the heat transfer coefficient (W m""^ °C"^), 

Nu is the Nusselt number, 

k is the thermal conductivity of the fluid (W m"^ °C), 

and 

D is the characteristic dimension of the system (m). 

The heat transfer coefficient, h^, is a value which has been averaged over 
the entire surface area of the system. The characteristic' dimension D 
must be defined for each different geometric shape one wishes to consider. 
Commonly, one might use the diameter of a sphere or the widest point of a 
leaf in the direction of the wind. The Nusselt modulus is a 
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nondimensional number used to scale laboratory results to other wind 

velocity and fluid properties. Kreith (1973, page 317) gives two physical 
^interpretations for the Nusselt number. It may be thought of "^as the 

ratio of the temperature gradient in the fluid immediately in contact with 
.the surface to a reference temperature gradient (Tg - T)/D" or the "ratio 

D/x where x is the fluid thickness of a hypothetical layer which, if 

completely stagnant, offers the same thermal resistance to the flow of 

heat -as the actual boundary layer/' 

In practice, another nondimensional number, the Reynolds number. Re, 
which is the ratio of inertial (pv^) to viscous forces (pV/D), is introduced 

to account for the scaling effects of fluid velocity, geometry, and fluid 
properties (see Kreith 1973 or Cowan 1977 for a discussion of Re). 




(12) 



where 

p is the fluid density (kg m~3), 
V is the fluid velocity (m s~^), 
D is the characteristic dimension (m), 

and 

p is the fluid viscosity (kg r?"^ s"^).* 

For a particular geometry, a functional relationship between the Nusselt 
number and the Reynolds number can then be calculated from laboratory 



* Sometimes the kinematic viscosity v .i^ s"^) is used in the 
Reynolds number where V = p /p • 
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measurements (h^ versus V). Usually, the result Is expressed as: 

Nu = a Re^ (13) 

^ where a and b are determined by regression. The parameters will be 
•^different for different shaped objects (Kreith 1973), then, using equation 
(12), the Nu can be calculated which in turn will yield the convection 
coefficient from equation (11). This procedure is the standard nethod 
used in engineering. In ecological studies, where one is concerned about 
the effects of wind speed and size and the geometry of the systen is 
constant, the convection coefficient is sometimes given without 
specifically indicating the coefficient of equation (13) as in the 
following exaniple. This approach is more direct and is fine to use as 
long as the student realizes what has been implied. 

Tibbals et al. (1964), froir. earlier experimental work by Gates and 
Benedict (1963), were able to use the following formula to evaluate 
convective exchaiige for a leaf; 

c = ^ (^) - V 

where 

C is the convective heat flux (W m"2 ^ positive for heat flow fror the 
leaf to the air) , 

kj is 9.1A (J m"2 oq-I s"1/2), 
V is wind velocity (m s"^), 

D is the characteristic dimension (widest point) of the leaf in the 
direction of the wind (m). 
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Tg is the surface temperature of the leaf (°C), 

and 

is the air temperature (*^C). 

After using this convection term C in the energy budget equation, 
Tibbals et al. (1964, page 538) concluded for similar energy environments 
**that • • • that the broad deciduous type leaf would be considerably 
warmer than the conifers . . and "... that the demand on transpiration 
is probably greatest per unit surface area for the broad~leaf plant than 
the conifers." (Also see Gates 1977.) 

The theory of convective transfer is more complex than presented 
here. Although some additional material will be included in the problems, 
the reader may wish to refer to Kreith (1973), Monteith (1973) or Campbell 
(1977). 

Evaporation 

Evaporation is the process of water changing from a liquid to a gas. 
For animals, water can be lost through respiration, through special 
glands, or through any part of the skin. Water loss is usually a small 
component of the heat balance for animals but may be large for animals with 
moist skins. It usually increases quickly at higher temperatures. Figure 6 
from Dawson and Templeton (1963), shows this effect for the collared 
lizard, Crotaphytus collaris* In general, the water loss, E, is given by 
equation (15). Again, this is in the form of a potential divided by a 
resistance (an Ohm's Law analogy). 

Co - (15) 
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Figure 6. Relation of evaporative water loss to ambient 
temperature in collared lizards weighing 25- 
35 g. Data represent minimal values at various 
temperatures for animals studied. (From Dawson, 
W.R., and J.R. Templeton. 1963. P. 231.) , 
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where 

E is water loss (kg s"-'), .< 
Cq - Cg is the water vapor concentration difference between the 
'surface (o) and the free atmosphere (a) (kg m"^), 

rg is the resistance to water vapor loss (s m""-^). 

The water vapor concentration at the surface of the water loss site may be 
saturated such as tiiat of a leaf or a lung of a bird but it may be less 
than saturated as in the skin of a lizard. The concentration in the air 
is equal to the saturated concentration (at air temperature T^ ) times the 
relative humidity. The resistance to water loss in the skin of a lizard 
will be controlled by the resistance to water movement in the skin and 
the boundary layer resistance. In a leaf, the total resistance to water 
loss is a combination of the stomatal resistance and the boundary layer 
resistance. But the water loss resistance from the moist skin of a frog 
was assumed to be equal to the boundary layer resistance only (Tracy 1976) . 
In general, the boundary layer resistance is a function of the diffusion 
coefficient of water vapor in air and the boundary layer thickness. This 
thickness in turn depends on the size of the object as well as the wind 
speed. 

To relate the water transport calculated in equation (15) to the 
energy balance, the mass flux E must be multiplied by the latent heat of 
evaporation L. L can be closely approximated as a function of 
temperature.* 

L(T) = 2.50 X 106^-^2.38 lO^ T, (16) 



* L is also a function of pressure. See Zemansky and Van 



Ness 1966« 
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where 



L(T) is the latent heat of evaporation (J kg"^), 



and 



T is the surface temperature at the site of evaporation (°C). 



-As an example of the energy flux, the water loss at = 40 °C for C 
collaris is 0.8 rag of H2O g""*'' hr"''' (Fig. 7). Correcting the units to 
kg s~^ for a 30 g lizard yields 6.67 x 10^^ kg s"^. Therefore, using 
equation (16), the energy loss for the aniraal is 1.6 x 10"^ W. In this 
example, we have not calculated E explicitly, mechanisticTally accounting 
for wind speed or water vapor concentration of the environment, but animal 
physiologists often report evaporation as a function of air or 
body temperature only. These factors become more important when water 
loss is a larger fraction of the total energy budget or when one is 
concerned with the water balance of the organism (Welsh and Tracy 1977). 
Later in the text, a functional relationship to include these factors is 
given for water loss fron a leaf. Generally, transpiration is a 
significant energy loss for plants which has led Monteith (1973) to make 
the distinction between "wet" and "dry" systems. 



In the above examples, we have mentioned several thermal properties 
such as thermal capacity and specific heat. Now we wish to define thep 
more carefully. 

Specific heat is defined as: 



THERMAL PROPERTIES OF MATERIALS 



c = 



AO , 

AT M 
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where 

c is the specific heat (J kg"^ °C"^), 
t AQ is the heat flux (J), 
•7 At is a 1° change temperature (°C), 

and 

.* 

M is a 1 kilogram system (kg). : 

We cannot measure the amount of heat a body can hold, so in a sense 
the term specific heat is misleading, but it is conunonly used. It is 
simply the amount of heat which must be added per unit mass of material 
for a 1° rise in temperature. What actually is changing is the internal 
energy of the system. Table 3 (after Monteith 1973) shows values of 
specific heat for different materials. See Stevenson (1979a) or Zemansky 
and Van Ness (1966) for a more formal definition. 

The heat capacity of the system can then be obtained by multiplying 
the mass of the system by specific heat: 

C = Mc. (18) 
For instance, if we consider 0.i2 kg of granite and the same weight 
of peat soil as two systems {see Table 3), the heat capacity of each is 
the specific heat c times the mass M or 96 J °C"^ for granite and 226 J 
°C~^ the peat. To calculate the energy necessary to raise each system 
10 ^C, we can employ equation (l8) in the form 

AQ = AT M c. . 

The answer is found to 960 J for the granite and 2,200 J for the peat. 
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Table 3. Thermal properties of natural materials. 
(From Monteith 1973) 



MATERIAL 


DENSITY 


SPECIFIC 
HEAT 


THERMAL 
CONDUCTIVITY 


THERMAL 
CAPACITY 


THERMAL 
DIFFUSIVITY 




P 

kg m 
X 10 


c 

J ks'-^'c"-" 

X 10"^ 


k 


c = pc 

V 

J .-3°c-^ 

X 10 


k 

pc 

2 -1 
m s 

X xu 


Granite 


2.60 


0.8 


4.61 


2.08 


2.22 


Quartz 


2.66 


0.8 


8.80 


2.13 


4.14 


Clay minerali^ 


2.65 


0.9 


2.92 


2.39 


1.22 


Ice 


0.9 


2.1 


2.30 


1.89 


1.22 


Old snow 


0.5 


2.1 


0.29 


1.05 


0.28 


New snow 


0.1 


2.1 


0.08 


0.21 


0.38 


Wet sand 
Dry sand 


1 A 
1*0 

1.4 


1 3 
0.8 


1.68 
0.17 


2.08 
1.12 


0.81 
0.15 


Wet marsh 
soil 

Peat soil 

Still water 

Still air' 


0.9 
0.3 
1.0 
0.001 


3.4 
1.8 
4.18 
1.0 


0.84 
0.06 
0.63 
0.07. 


3.06 
0.54 
4.18 
0.001 


0.27 
0.11 
0.15 
20.0 


Or Ecinic 
matter 


1.30 


1.92 


0.25 




1.00 


Fur 


0.98 










Mean body 
Fat 


0.98- 
1.05 


3.42 
1.88 


0.14- 
0.20 






Plant leaf 












Wood 


0.6 


1.3 


0.15 


.78 


0.19 


^ 
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An equally valid way of considering the heat capacity of a eysteir is 
to define the thermal capacity. It is the heat added per degree change 
per unit volume and thus is equal to the specific heat times tke density 
of the material of the system (see Table 3), 

AO 



where 



and 



Cy is thermal capacity (J oq-I m"*^), 

AQ is the heat change in the system (J), 
At is change in temperature (°C), 

V .is unit volume (m^). 



This is an especially useful concept if we wish to compare systems of 
equal volume. Let us reconsider our previous example of granite and peat 
except that now each system will have the same volume of 1,2 m^. Thus, 
they are not the same systems as before. The volumetric heat capacity for 
each is the product of the thermal capacity and the volume which for 
granite turns out to be 2.5 x 10^ J 0^-1 and for peat 0.65 x 10^ J °C"^. 
We have found an increase in the heat capacity for the second grafiite 
system and a decrease in the heat capacity in the second peat system. 
This result is expected if we consider the densities of each materia]. 

Thermal conductivity along with the notions of conductance and 
thermal resistance have already been discussed and defined in the section 
on conduction. The last parameter we wish to define is the thermal 
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diffusivity K which is the ratio of the thermal conductivity K divided by 
density p times the specific heat. 

k 



K = 

PC 



J m C s 



kg m-'j kg-' °C-' 



_ „2 ^-1 (20) 
= m s 

The thermal diffusity is important in nonsteady state conduction problems 
where thermal energy can be stored or transported as shown in the first 
example of heat flow in the soil. 

EXAMPLES OF THE HEAT ENERGY BUDGET 

To illustrate the heat transfer principles and their biological 
importance, we chose three systems: the soil, a leaf, and a lizard. 
Again, the reader is asked to think of the ecological importance of the 
system and how the thermal balance influences physiological processes as 
the examples are presented. 

Heat Flow in Soil 

In soil, heat flow has great biological importance. Not only does 
the energy balance control the soil temperature, but it greatly affects 
the moisture content. Both factors are important for photosynthesis in 
green plants. Soil temperature also influences nutrient uptake as we have 
seen in our example at the beginning of the module* Fungi, bacteria, 
insects, as well as many other invertebrates and vertebrates, commonly 
spend part or all of their life cycles in or surrounded by the^soil. Soil 
temperatures influence metabolic rates and may serve as a behavioral cue 
for daily and seasonal activity patterns. 
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The simplest representation of heat flow in the soil is a conduction 
process where the heat flux dQ per unit time dt is equal to the 
conductivity k times the difference in temperature dT over a unit path 
length dz. 



dQ ^ dT 
dt " "^dz* 



(21) 



The minus sign is used to indicate that the potential gradient is in the 

direction of decreasing temperature. 

Another result can be derived by considering Fig. 7. The First Law 

states that the change in internal energy (storage) is equal to the 
inflow minus the outflow. The change in internal energy will be equal to 

the nass (pdz) times the heat capacity (c), multiplied by the change in 

temperature with time (AT). 



dT 
dz 



- k 


fdT] 


A] 





AT = 



cp dz 

9T ^ k_ 
9t cp 



dT 




f y 

dT 




dz 

V, J 


U 


[dzj 





(22) 



Equation 22 is a form of the diffusion equation and is called the 
thermal diffusivity. This equation can be solved analytically if simple 
enough boundary conditions are assumed and the thermal diffusivity 
remains constant. In general, the diffusivity is a function of the water 
content of the soil and the soil temperature (here numerical methods must 
be used). We can, however, get a feeling for the heat transfer profile 
by noting the general boundary conditions. We know that at the soil 
surface under steady state conditions the temperature will cycle daily 
while at some depth 2^ below the surface the temperature will not change. 
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T' 0* K {dT/d7 ) 




Outflow k IdT/dr ^ 



Figure 7. Schematic view of a unit of soil, 
showing heat inflow from above 
and heat outflow below, leading 
to derivation of Eq. 22. (From 
Lowry, W.P. 1969. P. 51) 
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Figures 8 and 9 from Lowry (1969) give different representations of how 
the soil might change with time and depth. Simpson (1977), Carslaw and 
Jaeger (1959), and Van Wijk (1966) present a more detailed disQission of 
soil heat flow. 

A Leaf 

Our second example is the energy balance of a leaf. Gates (1968) 
chose to model the heat balance of a leaf because it has a distinct 
physical geometry and because it is the unit of photosynthesis. The 
governing heat flux equation for the steady state condition was given in 
the module on the First Law of Thermodynamics and is repeated here: 

Qa = Qe + C + LE, (23) 

where 

Qa is the radiant energy absorbed by the leaf (W m"2)^ 
Qe is the energy emitted by the leaf (W m"^), 
C is the convective flux (W m"2)j 

and 

LE is the evaporative flux (W m""^). 
With empirical measurements and some assumptions (see Gates 1977), 
equation (23) is expanded as: 



- ca(T^ 273)' + 



0 - 

V 



C (T ) - r.h. C (T ) 

where 

a is the Stefan-Boltzmann;: constant (5.67 x lO^^ y ^-2 k~^), 
e is the emissivity of the leaf, 
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Ttmperoture 



Figure 8. Generalized soil-air temperature 
profiles (tautochrones) near the 
soil surface, for four-hour inter- 
vals during a diurnal period. 
(From Lowry, W.P. 1969. P.37J 




Figure 9. Generalized diurnal patterns of 
isotherms near the soil surface 
on coordinates of time and 
distance. (From Lowry, W.P. 1969. 
P. 37.) 
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is the leaf temperature (^c), 
V is the wind speed (m s~^), 
ki is an experimentally determined coefficient (9.14 J m"2 - 
s-0.5 oc-1), 

k2 is an experimentally determined coefficient (183 s""0-5 m"^), 
Tg is the air temperature (°C), 
D is the leaf dimension in the direction of the wind (m), 
r.h. is- the relative humidity, 

W is the leaf dimension transverse to the wind (m), 
Cq(Tj^) is the concentration of water vapor saturation at leaf temperature 
(kg 

C (T ) is the concentration of water vapor in free air saturated at air 
a a 

temperature T (kg m 
a 

and 

r^ is the internal diffusion resistance of the leaf (s m"^). 

Inspection of equation (2A) shows that the energy balance of the 
leaf is very dependent on leaf temperature which in turn affects 
convection and transpiration. To analyze this multivariate problem, 
Gates graphed rhree variables while holding the rest constant. Figures 
10 and 11 illustrate thermodynamic conditions that would be common in 
tropical habitats at midday sunny (Q^ = 1.2 cal cm"2 ^in) and cloudy 
(Qa ~ cal cm"2 min) radiation levels. Note the differences in 
transpiration rates and leaf temperatures as a function of leaf 
resistance. Gates (1968) also considered the effects of leaf size. 
Figure 12 clearly shows that small leaves, by reducing leaf temperature 
and transpiration, would be advantageous in desert regions. 
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LEAF TEMPERATURE. •€ 



Figure 10. Computed transpiration -rate versus leaf 
temperature as a function of wind speed 
and internal diffusive resistance for the 
conditions indicated . 
(From Gates, D. M. 1968. P. 216.) 
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Figure 11. Computed transpiration rate versus leaf 
temperature as a function of wind speed 
and internal diffusive resistance for t.h 
conditions indicated. \ 
(From Gates, D- M. 1968. P. ^17.) 
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LEAF TEMPERATURE, 



Figure 12. Computed transpiration rate versus leaf 
temperature for various leaf dimensions 
and internal diffusive resistances for 
the condition:5 indicated. 
(From Gates, D. M. 1968. P. 229.) 
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A Lizard 



As a final example, we will review the work of Bartlett and Gates 

.< 

(1967) who computed the heat energy budget for the Western fence lizard, 
.Scelo^orus occidentalism on a tree trunk. They hypothesized that 5. 
oocidentalis oriented itself on the tree trunk so as to maintain its bod 
temperature relatively high and constant. The daily temperature 
variation of the tree surfaces supported their idea (Fig. 13). The heat 
energy equation for equilibrium is 

Qa + M - Qe - LE - C - G = 0 (25) 

where 

Qg is energy absorbed (W m"^)^ 

M is metabolism (W m"^), 
Qg is radiation emitted (W m"'- ) , 
LE is evaporation (W m 

C is convection (W it ^ ) , 

and 

G is conduction (W Tn"^)^ 

Equation (25) may be rewritten: 

+ M - C a Ag (273 + Tg)^ - LE 
-h. A, (T3 -T,)-kAe^ = 0 (26) 

where 



e Is emissivity to longwave radiation, 
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Figure 13. Tree surface temperature in °C 

on Chew's Ridge at 10:00, 1.^:00, 
and 17:00 Pacific Standard Time 
June 21 as a function cf direction 
from North. 

(From Bartlett, P. N. and D. M. 
Gates. 1967. P. 320.) 
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a is Stefan-Boltzmann constant (5.67 x lO"^ W in'**^K"'^), 
Ag is the percent of the total surface area of the animal not in 
contact with substratum, ^ 
Tg is surface temperature of the lizard (°C), 
is convection coefficient (W m"^ oc""l), 
is air temperature (°C), 
k is thermal conductivity of substratum (W m"^ °C"^), 

is the percent area of the animal total surface area in contact 
with the substratum, 

^ is temperature gradient with the substrate (°C rrT-^). 

The terms of equation (26) were either experimentally determined or 
taken from the literature. It was then solved for a variety of 
situations to obtain a range values for the animal's energy balance: 

1. A maximum heat gain: If the lizard's surface temperature Tg 
equaled 29.3 °C and the lizard was in close contact with the tree. 

2. A minimum heat gain: If Tg = 39.9 and the lizard was not in 
contact with the tree. 

3. Maximum heat loss of Tg= 38.9 °C, wind speed = 2.23 m s"^ (5 
mph) and the animal was oriented at 90° to the direction of the wind. 

4. Minimum heat loss if Tg - 29.3 ^C, and wind speed was 0.1 m s""^. 

5. Intermediate heat loss if Tg = 34.1 and wind speed = 0.45 m 

s-1. 

Figure 14 gives sample results. 

Finally, they calculated the energy losses and gains every hour at 
45° increments around the tree. Using the criterion that the energy 
losses must fall between the energy gains (Fig. 14), they could predict 
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PACIFIC tT»NO»BO TimC 

Figure 14. Max and min energy gains (concave downward) and 
Tnaximum, miniTnuin, and estimated energy losses 
(concave upward) of a lizard on a tree trunk on 
Chew's Ridge for selected times and positions on 
June 21, with times when lizards could be ex- 
pected to be found at the selected positions. 
Areas hatched represent times when the estimated 
energy loss falls between maximum and minimum 
energy gains. (From Bartlett, P. N. and D. M. 
Gates. 1967. P. 321.) 




X.MtUAL MIIRVIt »«*n»«RI 



Figure 15. Estimated hourly positions on p tree trunk on 
Chew's Ridge for June 21 where Sceloporus 
occidentalis could maintain its characteristic 
body temperature (arrows) and actual observed' 
positions of lizards on June 21, 196A (crosses) 
(From Bartlett, P. N. and D. M. Gates. 1967. 
P. 321 .) 
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the position of the lizard on the tree trunk as a function of time of 
day* The predicted and observed positions of the lizards are compared in 
Fig. 15. .< 

This module has attempted to develop the aspects of heat transfer 
'that are important to biology* Emphasis has been placed on the physical 
processes, radiation, evaporation, conduction and convection, that 
influence the physiological, behavioral and ecological activities of all 
organisms. The discussion in conjunction with the modules on the First 
Law of Thermodynamics presents the basic physical laws that are needed to 
describe the heat energy balance of ecological systems. 
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PROBLEMS 



lat To measure the oxygen consumption of resting homeotherms as a 
^function of the thermal environment, animals are commonly placed in an 
environmental chamber where the air temperature can be changed and 
•monitorod. The results of the experiment are given as a plot of O2 
consumption versus air temperature. Imagine that for some reason the 
wall temperature is not equal to the air temperature. Plot the radiation 
emitted by the wall per unit surface area, Qew> a function of wall 
temperature using Stef an-Boltzir.ann ' s Law (assume emissivity equals 0.96). 

lb. If the emissivity of an environmental chamber with plastic 
walls is 0. 6 for longwave radiation what is its absorptivity to longwave 
radiation? Imagine that the air is much cooler than the animal's 
surface. The radiation leaving the animal, incident on the walls, could 
heat the walls above the air temperature and increase the radiation 
incidence on the animal. Name two ways to avoid this problein. A chamber 
with metal walls will have a small emissivity ('X'.IS) to longwave 
radiation. What will be its reflectivity? What will happen to the 
radiation leaving the animal in this case? How cm this effect be 
minimized? 

2a. It is also possible to represent Planck's Law as a function of 
frequency f rather than wavelength X; 



27^h 



r3 



exp (— ) -1 



-1 



Equation (1) of the text gives the relationship between A and f . Show how 
Ef and E^ are related. What is the difference in dimensions between and 
E^? The concept of wave number is often introduced to plot E^ and F.-)^ of the 
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same graph. Wave number n = X"'^. How are f and n related? Plot E^^ and 
Ef for T = 280 K as a function of X. Repeat for T «= 6,000 K but show both 
the X and n scales. How does the relationship of bandwidth/ w^^e number 
and bandwidth/wavelength change over the spectrum? 

2b. Planck also found that the energy in a photon, e, is 
proportional to its frequency (e = hf). What does this imply for a E^ 
versus Ef representation? At what wavelength does the sun radiate the 
mos ,^y? Assume its surface temperature is 6,000 K. How does this 

value compare with that in Fig. 4? 

3a. Using equation (8), plot the heat of vaporization as a function 
of temperature. Schmidt-Nielsen (1975, page 313) says that 580 cal g"^ 
is a commonly used value in physiology. Is this a good approximation? 
Evaluate the error over the range 0 to 40 °C. 

b. Schmidt-Nielsen (1975, page 37) gives the lung volume for 
mammals as a function of weight as 

V = 0.0567 m1-02. 

V is lung volume (liters) and M is body mass (kg). 

If respiration rate (breaths rain"-) is a function of air temperature 
(for a moose (Belovsky 1978)), 

Br = 7.72 exp (.07 T^), 
calculate the respiration loss in watts as a function of air temperature 
(Ta = -20 to 40 °C). Assume M = 358 kg, that the exhaled air is at air 
temperature, the relative humidity is 50 percent and L is constant at 
2.43 X 10^ J kg"^. Repeat the calculation for the exhaled air at 

T + T, ! 
a b . 

2 ■ ' 
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Hiot: The energy lost is equal to the net mass flow (the mass of 
water exhaled minus the mass of water inhaled per unit time multiplied by 
the latent heat of water). Convert the volume exchange to a id«ss 
'exchange using the fact that concentration (mass/volume) of water is the 
"^product of the molecular weight, the vapor pressure of water at that 
temperature and the relative humidity divided by the product of the gas 
constant and the absolute temperature of the sample. 

4. The rate of heat transfer by conduction depends on the 
potential (temperature difference), the area of contact and the 
conductivity of the material. Estimate the heat transfer between an 
animal and the ground for all combinations. 

Surface temperature = 20 ^ 
Ground temperatures = 2, 14, 27 **C 
Area of contact = 0.054 m"^ 

Thermal conductivities = 0.02, 0.08, 0.17 W m~^ ^C"^ 
Thickness of the material layer limiting heat flow = 0.005 m 

5a. In the diagram, heat flows through two materials from the left 
boundary (temperature Tj ) to the right boundary (temperature T2). Assure 
steady state conditions. 



49 



45 



0) 

u 
o 

CO 
u 

0) 

a 
B 
0) 
H 




How does the temperature drop at the boundary depend on thermal 

ki M 

conductivity and the thickness of the material? Let , = a and * ~ b; 

TB-T2 2 
find Tg. Let ^ and find Tg. Plot" ^ _ If a = 1.2, b = 0*2, Tj 

= 40 °C and 10 °C, in which material is the temperature gradient 

larger? Which material has the larger temperature drop? What is Tg? 

b. In the conduction section of the text the heat flow between an 
alligator and the substrate was calculated. It was assumed that the heat 
flow into the rock was not important. Discuss this assumption with 
regard to part "a" of this problem. 

c. Imagine that the alligator's skin temperature is initially 
equal to the body temperature. How will the temperature profile change 
in the soil and in the animal as the heat exchange .approaches a constant 
value? 

6. Kreith (1973) gives the- following relationship for Nu and Re of 
cylinders: r ' 

Nu = a Ke^ 
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Re a b 

O.A-A 0.891 0.330 

4-AO 0.821 0.385 

AO-A,000 0.615 0.A66 

A,000-AO,000 0.17A 0.618 

AO,000-AOO,000 0.0239 0,805 



Assume that wind velocity is 2 m s"^ and 0.1 m s"^, the diameter is •A6 m 
and .002 m. Compute h^ for the four combinations of diameter and wind 
velocity. Take V = 1.A2 x 10"^ m^ s"^ and k = 2.50 x 10^2 w a"-^ K^"^ at an 
air temperature of 10 °C. 

7a. Calculate the heat transfer coefficient for a flat plate using 
Nu = 0.60 Re^*^ (Monteith 1973, page 22A; Kreith 1973, page 3A1). Let 
air temperature be 20 ^C, so V = 1.51 x lO^^ ^^2 g-l ^ 2.53 x 10^2 ^ 

ijj-1 K"l, How does this compare with the convection coefficient for the 
leaf? This is accurate for laminar flow Re < 2 x 10^. 

b. Show the region Re < 2 x 10^ on a graph of wind velocity V 
versus the characteristic length D. When Re < 2 x 10^, Monteith suggests 
0.032 Re^'^. How likely is this to occur naturally? Would a leaf be 
rigid in these wind speeds? 

8. In the text, we discussed the physical process of free 
convection. To determine the importance of free convection, one must 
compute the ratio of the Grashof number to the square of the Reynolds 
number. If the ratio is approximately 1, then free convection' cannot be 
ignored. The Grashof number is • ' . 
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. agD^d - T) 
or.. 1 

V 

g ^ gravitational acceleration (9.8 m s""^)^ ^ 

a = coefficient of thermal expansion for fluid (a = 1/273 for air), 

Tg = surface temperature (°C), 

T = fluid temperature (^C) 

V = kinematic viscosity (m s ). 
<;r 

Plot ^ on a graph of V versus D, Let - T equal 10, 20 and 30. 
Re^ ^ 

Below the line Gr = Re^, free convection loss cannot be ignored. Is free 
convection important for^ leaves in the natural environment? 

9a. The energy balance for fruits is important in both cultivated 
and wild plants because it affects when and at what rate the crop ripens. 
Grapes for instance must be picked at just the right time to insure the 
>QSl sugar content for winemaking. Mcintosh apples need cold nights and 
warm days to be of the best quality. When weather conditions are 
unfavorable, fruit can become overripe or rotten before it is harvested. 
Wild plant fruits are an energy investment by the plant to attract seed 
dispersers. This is especially true of birds and bats in the tropics. 
The color of the fruit is often an important signal to a potential 
disperser while the size and spacing of the fruits will influence the 
animal's harvesting rates. These same characteristics are important to 
the fruit's thermal balance. 

An empirical approach based on the "heat unit" or "growing degree 
day" concept has been developed to predict the harvest times. When the 
fruits have been exposed to x degree days the theory says they will be 
ready to pick. Degree days are computed as follows: • 



i. 
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DD = Z 



■^amax ~" '^amin 



2 



where DD is degree days, T^jjj^j^ is the maximum air temperature during the 
day, Tgjjj^^ is the minimum air temperature for the day and T is^a threshold 
^temperature specific for each species and which can be thought of as the 
temperature below which development is stopped. (The summation is over 
only positive values of the term). From the data below calculate the ri- 
pening degree days if is 16 **C. 

Day 
'^amax°^ 
'^arain°^ 



1 


2 


3 


4 5 


6 


7 


8 


9 


10 


20 


25 


21 


21 18 


20 


25 


27 


26 


23 


16 


18 


U 


13 10 


12 


15 


18 


18 


16 


are 


there 


with 


this approach? 













b. A more mechanistic viewpoint should allow one to calculate the 
temperature of the fruit. As an example let us consider a spherical 
model 1, 4 and 10 cm in diameter. What is the heat energy balance 
for the fruit. Which terms are small and can be ignored? What terms are 
you unsure about? As a first approximation, let us start with the 
following model: 
Qa = Qe + C, 

where 

is absorbed radiation (W m""^), 
Qq is reradiation (W m"^), 

and 

C is convection (W m^^)^ 
Under steady state conditions the temperature of the fruit core, T|:, 
will equal the temperature of the surface. Therefore, = e a Tf^ and 
C = h^ (Tf - Tg^), For forced convection about a sphere without the 
increased turbulence of the outdoor environment (Nobel 1975; Kowalski and 
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Mitchell 1976; Monteith (1973, page 224) give two formulae for the 
relationship between the Nusselt and the Reynolds number. 

Range of Re Nu 

0-300 2 + 0.54 ReO-5 

50 - 1 X 105 0.34 Re^-^ 

Calculate Nu by both formulae for Re = 50, 100 and 300. If D = 1 cm, 
what is the minimum wind speed for which the second relationship is 
valid? Which formula yields a lower heat transfer coefficient? What will 
the effect on the energy balance of the fruit be? Find h^ as a function 
of D (0.01, 0.04 and 0.1 m) and V (0.1, 1.0 and 10.0 m s""M if Nu = 0.34 
Re^*^. Assume v = 1.42 x 10"^ m^ s"^ for the entire problem. 

For a sphere: = [AjS + A2S + A3r +(S+s)] + ^a^a^^a + r^l 

where 

ag = the absorptivity to shortwave radiation (assume .7), 

Aj = the percent of the total surface hit by direct beam radiation 

2 2 

weighted by the angle of incidence (equals Tir /47Tr = 1/4) 
< 

S = direct beam solar radiation (W m"^), 
^2 = Ag = A^ = the percent of the total surface hit by the 
corresponding radiation sources (assume 0.5), 

s - diffuse shortwave radiation (assume = 0.1 of total shortwave 
radiation Rp , W m"^) 

r = reflectivity of the ground (assume 0.15), 

- absorptivity to longwave radiation (assume 1.0), 
R3 = longwave radiation from the atmosphere (W m^^), 
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and 

Rg = longwave^ radiation from the ground (W Tn"'^). 

From Morhardt and Gates (1974, page 20, Fig. 2c), I have tTaken R-,, 
-R3 and Rg for a sunny day and converted the units to W in"'^ tabulated 
jbelow: 

TABLE FOR PROBLEM 9b 



Hour of 
the Day 


R 

P 


R 

a 

(W m"') 


R 

g 

(W ni"^) 


(°C) 


7 


35 


223 


335 


7.5 


8 


593 


223 


3A8 


10.0 


9 


837 


230 


369 


12.0 


10 


loi: 


237 


A19 


15.0 


11 


1116 


lUU 


A61 


19.0 


12 


1186 


251 


502 


24.0 


13 


1116 


253 


A81 


27.0 




1012 


251 


A 61 


25.0 


15 


837 


lUU 


419 


22,0 


16 


593 


237 


391 


16.0 1 


17 


35 


237 


363 


12.5 


18 


0 


230 


3A9 


10.0 



Calculate values for each hour of the day. What part of this is 

shorturave radiation? 

Assuming that the wind speed is 1.0 m/s and that the diameter of the 

fruit is 0.04 m and using the convection coefficient you have derived, 

calculate the fruit's equilibrium temperature for each hour of the day. 

Using and T values for 7 and ? o'clock calculate the equilibrium 
a a 

fruit temperature for all co,nbinations of V = 0.1, 1.0 and 10.0 m/s and D 
= 0.01, 0.04, 0.10 m. Size and wind speed affect the fruit temperature 
in opposite ways at the two times. Why is this? What are the rptoblems 
with this approach? 
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10. Discuss the following group of papers concerning heat transfe 
in biology: 

'Strunk, T. H. 1971. Heat loss from a Newtonian animal. J. Theoret. 
Biol. 33:35-61. 

Kleiber, M. 1972. A new Newton's Law of cooling? Science 
178:1283-1285. 

Tracy, C. R. 1972. Newton's Law: its applicability for expressing he 
losses from hpmeotherms. Bioscience 22:656-659. 

Technical Comments: Perspectives of linear heat transfer: T. H. Strunk 
' C. R. Tracy, M. Kleiber- Science 181:184-186. 

Linearized heat transfer relati'>. s i?i biology. G. S. Bakken and D. M. 
Gates; T. H. Strjnk; M. Keiber. ^Science 183:976-978. 

Bakken, G. S., and D. M. Gates. i974. Notes on "Heat loss from a 
Newtonian animal/" J. Theoret. Bid. 45:283-292. 
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PROBLEM SOLUTIONS 

1. 




la. Metabolic heat will increase the wall temperature. The smaller 
the chamber, the greater this effect will be. The influ/nce will change 
as the square of the distance. 

b. The absorptivity is 0.96. Two ways to minimize wall heating 
are to make the chamber large relative to the animal or immerse the 
chamber in water. If the emissivity for the metal walled chamber is 0-15 
absorptivity will be 0.15. Because the walls will not transmit longwave 
radiation, 85% will be reflected back to th^ animal and the other sides 
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of the chamber, again making the chamber feel wanner to the animal than 
the air or wall temperatures would indicate. Changing the inside surface 
: chaidcceristics of the chamber cr increasing the chamber size relative to 
.the animal will reduce the importance of reflected radiation. {JSee 
Porter 1969 and Morhardt and Gates 1974 for discussion.) 

2a. Ef = X"2 E^/cL. The dimensions of Ef are M L~"2 £^ ^ 

1^3. They differ by a factor of L^^ T^^, If X f = c^ and n = X"^ 
then n = f/cL. (For plots see data below and the following figures.) 
Because wavelength is inversely proportional to wave number, equal 
bandwidths on a wavelength plot will cause unequal bandwidths on a wave 

Table for Problem 2. 



T 


A 


n 


I 


A 


z 






X lU j 


V.5 X 10 j 


(J m s } 


(J m ) 


280 


U,0 


0.250 


0.750 


0.10 X 10^ 


0.05 3: 10"^^ 




5.0 


0.200 


0.600 


0.41 


0.35 




6.0 


6.167 


0.500 


0.92 


1.10 




7.0 


0.143 


0.429 


1.45 


2.37 




8.0 


0.125 


0.375 


1.86 


3.98 




9.0 


0.111 


0.333 


2.11 


5.71 




10.0 


0.100 


0.300 


2.21 


7.38 




11.0 


0.091 


0.273 


2.20 


8.87 




12.0 


0.083 


0.250 


2.11 


10.31 




13.0 


0.077 


0.231 


1.98 


11.14 




15.0 


0.067 


0.200 


1.66. 


12.45 




17.0 


0.659 


0.176 


1.35 


13.00 




20.0 


0.050 


0.150 


0.97 


12.95 




30.0 


0.033 


0.100 


0.34 


10.17 




AO.O 


0.025 


0.0/5 


0.14: 


7.46 




60.0 


0.017 


0.05? 


0.04 


4.25 




100.0 


0.010 


0.030 




1.85 




150.0 


0.007 


0.020 




0.90 


6000 


0.15 


6.66 


20.00 


0.06 X 10^^ 


X 10"" 




0.20 


5.00 


15.00 


0.73 


0.10 




0.25 


A. 00 


12.00 


2.63 


0.55 




0.30 


3.33 


10.00 


5.22 


1.57 




0.35 


2.86 


8.57 


7.57 


3.09 




0.40 


2.50 


7.50 


9.15 


4.83 




0.^5 


2.22 


6.67 


9.91 


6.69 




0.50 


2.00 


6.00 


10.00 


8.33 




0.55 


1.82 


5.46 


9»65 


9.73 




0.60 


1.67 


5.00 


?.03 


10.83 




0.65 


1.5A 


4.62 


£.28 


11.66 




0.70 


1.43 


4.29 


7 .50 


12.24 




0.80 


1.25 


3.75 


C»01 


12. 8i 




0.90 


1.11 


3.33 


4 74 


12.82 




1.00 


1.00 


3.00 


3.75 


12.48 




1.25 


.80 


2.40 


2.1' 


11.00 




1.50 


.67 


2.00 


1.25 


9.37 




2.00 


.50 


1.50 


.51 


6.73 




3.00 


.33 


1.00 


.13 


3.77 




5.C0 




.60 


.02 


1.62 




10.00 


.10 


.30 


SU— 


.46 
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number plot. For equal wavelength bandwidths larger wave numbers imply 
larger bandwidths. 

t 

b. A frequency or wave number plot has the property that equal 
areas represent equal amounts of energy which is not true of a wavelength 
plot. The maximum energy is at A =^ 0.85 x 10^6 ^ which is in the near 
infrared portion of the spectrum, while Fig. 4 shows \^^^ ~ 0.48 x 10"^ 

3a. Equation (16) is 

L(T) - 2.50 X 10* - 2.38 x lO^T 



2.5 



2.i* 



2.3 



30 



>.i3 



580 cal g"l = 2.43 x 10^ J kg^l . 
Body temperature for most homeotherras is above 36 ^c. Expelled air will 
be below body temperature so this seems like a good estimate. If the air 
was really at 0 ^C, the error would be: 



2.50 - 2.43 



2.9 percent, 



2.43 

At 45 °C, the error would be: 



2.43 - 2.40 
2.43 



= 1.2 percent, 
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b. 



For the moose 



Vl = 0.0567 (358)^-02 = 22.83 



liters 
breath 



U.iters of air per min are: 



c = 176.3 exp (.07 T3), 

Ry = 2.94 exp (.07 T3) in3 10^35-1. 

Assume that the air is saturated at T^^ = 38 ^C* We need to calculate the 

kg of water vapor in a liter of air (C^, kg m"3) exhaled at (T^ + 13)72 
^ me(T) 
~ RT 



e(T) is saturated vapor pressure of water at temperature (T) (mbars). 



where 




and 



18 g mol" e (T) 



C = 



8.31 kg m^s"^mor-^K"-^ T 



e 



.217 e (T + T )/2 



(T, + T )/2 
D a 



The concentration of water vapor inhaled is 
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LE 



2.43 J kg-1 (Cg - Ci) • 



2.43 r.h. 



.217et(T + T,)/2] .217e(T )r.h/ 
a b - a 



(T. + T )/2 
b a 



X 2.94 exp (.07 T ). 

a 



Table of evaporative water loss as a function of expelled teoperature. 



















X 10-' B-* 


u 




T +T, 


















T 

a 


a b 
2 


C 

_e 


C 
e 


fi 




e i 


B 

_ V 


LE 


LE' 


-20 


9 


.0090 


.001 


.0005 


.0085 


.0005 


.73 


15 


0.90 


-10 


14 


.0120 


.002A 


.0012 


.0108 


.0012 


1.45 


38 


4.2 


0 


19 


.0164 


.0068 


.0026 


.0140 


.0026 


2.94 


100 


17.1 


10 


24 


.0218 


.0096 


.00/. 8 


.0170 


.0068 


5.92 


264 


69.1 


20 


29 


.0290 


.0172 


.0086 


.0204 


.0086 


11.92 


591 


269 


30 


34 


.0375 


.0302 


.0151 


.0224 


.0151 


24.01 


1307 


881 


40 


39 


.0685 


.0510 


.0255 


.0270 


.0255 


48.34 


3172 


2996 



4, The conduction transfer during steady state can be calculated 
according to the following formula: 

kA (T^ - T ) 



G = 



c b g 
Ax 



Table of G values (W) 



T^C^kCW m-^ °C-^) .02 ^8 ^7 

2 3.89 15.55 33.05 

14 1.30 5.18 11.02 

27 -1.51 -6.05 • -12.85 
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DO 



5. By definition 

k k 

^From the First Law of Thermodynamics, we know that 

k k 

5 ^1 = "^2- Therefore, — (T^-Tg) = ^ (T^-T^) 



Solving for Tg, we have 





+ T 


M 






L 

i iJ 


2 


,"2j 







T = T 
B ^1 



let k,^ = ak2 and' L^= bL^ , then 



+ T, 



bL^k^ 



L^ak2+ bL^k^ 



T„ = 



L^k, 



B L^k^ 



^1^ 



a + b 



^2 ^ 



a + b 



Now let = X 
b 



xfl ('^l^ ^2^ 



Since a > 1, > k^ so material 1 has the steeper temperature gradient, 
a 1.2 

Because x = ~ = - = 6.0, the total temperature drop is greater in 



material 2. = 35.7 

B 
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b. A more realistic description of the heat transfer would allow 
for a slight temperature drop ih the soil. Then, in part a , « T and 
^2 = Tb- 

c. The skin will warm from its initial temperature of 22 oc 
towards ground temperature. Because heat is leaving the soil there will 
be a thin layer that cools down from 30 "C. Eventually an equilibrium 
will be reached and the temperature profile will be constant. 

6. 



Re 


V 








O.I 


2.0 


.46 


3239 


64789 


1002 


' 14 


282 



Nu 



Nu = a.Re 











0.1 


2.0 


.46 


26.59 


178.51 


.002 


2.27 


8.52 



Nu k 
D 



D ■ 



0.1 



2.0 



.46 . i 

r 

.002 '■■ 



1.45 

28. 38 



9.70 
106.5 



7. 



Nu = 0.60 Re°-^ , h = 



C 

10.5 



Nu k 
D 

0-60 , fv , 

7^0075 ^ Nd 



0.5 



= 3.98 

c 



V 

Di 



0.5 
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This is approximately 2.3 times the value of 9.14 J °C"^ s"^/2 
that Gates uses. Most wind speeds in natural environments are less than 
3 ra s"^ and most leaves are less than 0.1 m along their largesJ^ 
* dimension, thus, Re < 20,000, Leaves usually will not flutter-at wind 
^speeds less than 3 m s"^, {See figure for Problem 7 on following page.) 
See Grace (1977) for additional information. 



Figure for Problem 7 



Figure for Problem 8 



Curves for which the Reynolds nunber squared cqujls :hv Crasnof 
nutabcr . 



1 

t, > 



□ s *) 



.1 .2 .3 

LMawcti-r of lt,if in cx-iers 



.5 



.A 



.3 




0.1 0. 
Dlaneifr (-} 



8. 



Re2 = Gr 



Re = Gr 



V = 



ag D(Tg - T) 



V = 



.1895/(T - T-)D 
s 



61 



Velocity (m 5 ) as a function of - T and D 
Ts - T CO 



D(m) 


5_ 


10 


20 


.3 


.232 


.328 


.568 


.2 


.189 


.268 


.464 


.1 


.134 


.189 


.328 


.05 


.095 


.134 


.232 


.01 


.042 


.060 


.104 



From Figures 10, 11, and 12 in the text it is clear that leaf temperature 
is usually not more than 10 °C above air temperature (although this is a 
function of leaf size in Fig. 12). Leaf size is usually less than .05 m 
and wind speed is below 3 m s~^. There still remains a set of conditions 
where free convection may be important. For additional information see 
Tibbals et al. (1964) and Grace (1977). 

9a. The degree devs DD = 30. This approach assumes: 1) that the 
high and low air temperatures will represent the average thermal 
conditions that '.ihe fruit experiences; 2) that development is a linear 
function of the high and low average; and 3) that there is no problen? 
with the fruit overheating. There are undoubtedly other problems also. 

b. For Re = 50 and D = 0.01 m, 1.42 x 10"^ 

DV . . „ Re V 50 x 1.42 x lO'^ , 
Re = — So V - — « — = .710 m/s 

A comparison of the Nusselt numbers at the same Reynolds number (in 

the range 50-300) shows that the Nusselt number computed by the first 

formula is larger. Since th^ heat transfer coefficient is proportional 
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to the Nusselt number, the choice of the second formula will give a lower 
convection coefficient. Choosing the second formula will make the 
fruit's temperature depart from air temperature to a greater eftent than 
the first formula. 

0.34 (^)°-^ k 

Re=5V 3„ Nu =0.34(f)°-' and = ^ • 

Therefore, if 

k = 0.025 W m^"^ k"""^, h = 6 . 8 V° * ^ * ^ * 

c 

h for a sphere 
c 



VCm/s) 



DCm) 


0.1 


1.0 


10.0 


0.01 


10.8 


42.9 


170.8 


0.04 


6.2 


24.6 


98.1 


0.1 


4.3' 


17.1 


68.0 



Qa values are in the table below. Because ag, A]^, A2, A3, r and S/s 
are assumed constant, the absorbed shortwave radiation is equal to 0,245 



7 


8 


9 


10 


11 


12 


13 


li 


15 


16 


17 


^3 


238 




505 


576 


626 


667 


6:. 2 


60i 


536 


^59 


30S 


290 


7.5 


10.0 


12.0 


15.0 


19.0 


2^.0 


27 


25 


22 


16 


12.5 


10 


6.0 


12.9 


17.0 


21.9 


26.8 


32.1 


33.7 


30.9 


26.2 


18.8 


10.5 


8.1 



For 7 o'clock 



Tf values (°C) 

V(m/s) 



D(m) 


0, 


,1 


1.0 


10.0 


0.01 


4, 


.6 


6.6 


7.2 


0.04 


3, 


.3 


6.0 


7.1 


0.1 


2, 


,4 


5.4 


6.9 
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For .2 o'clock 

Tf values (^C) 

V(m/s) 



D(m) 


0.1 


1.0 


10.0 


0.01 


38.7 


29.1 


25.4 


0.04 


43.9 


32.1 


26.4 


0.1 


47.3 


34.8 


27.4 



The Tc values given can he more accurately computed for the 
conditions given but the estimates of Rp, , Rg and do not warrant 
it. If Lhe absorbed radiation is below that found in a blackbody cavity 
at air temperature Tg, then decreasing the convection coefficient 
(decrease V or increase D) will increase the difference between and Tf 
making the fruit colder. Conversely if is greater than a blackbody 
level increasing the convection coefficient will bring the fruit closer 
to air temperature. The problems with this approach are several. 
Usually the radiation fluxes are much more complex than we have 
considered because of surrounding vegetation or other fruit. Clumping of 
fruits will also decrease the convection coefficient. The microclimate 
on opposite sides of a plant is often different and difficult to specify 
accurately in small places. 

The fruits grow and change color during the season. Smart and 
Sinclair (1976) have investigated some of these problems as well as the 
temperature differences within one fruit and nonsteady state conditions. 

10. No solution provided. 



69 



64 



LITERATURE CITED 

Bartholomew, G. A. 1977. Body temperature and energy metabolism. Pages 
364-448 in M. S. Gordon, ed. Animal physiology: Principles and 
adaptation. Macmillan Publ. Co., New York. ^ 

1 

Bartlett, P. N., and D. M. Gates. 1967. 'ihe energy budget of a lizard on 
^ a tree trunk. Ecology 48:315-322. 

•"Belovsliy, G. E. 1978. Life history strategies of the moose. Ph.D. 
D.: sseruation. Harvard Univ., Cambridge. 

Boylen, L. W. , and T. D. Brock. 1973. Effects of thermal additions from 
tne Yellowstone geyser basin on the benthic algae of the Firehole 
River. Ecology 54:1282-1291. 

Byers, H. R. 1974. General meteorology. McGraw-Hill Book Co., New York. 
461 pp. 

Campbell, G. 1977. An introduction to environmental biophysics. 
Springer-Verlag, New York. 159 pp. 

Carlson, D. M. , and J. N. Gentry. 1973. Effects of shading on the 

migratory behavior of the Florida harvester ant Pogonomymex badius. 
Ecology 54:452-453. 

Carslaw, H. S., and J. G. Jaeger. 1959. Conduction of heat in solids. 
Clarendon Press, Oxford. 510 pp. 

Chapin, F. S. 1974. Morphological and physiological mechanisms of 

temperature compensation in phosphate absorption along a latitudinal 
gradient. Ecology 55:1180-1198. 

Collins, J. C, T. C. Pilkington, and K. Schmidt-Nielsen. 1971. A model 
of respiratory heat flow in a small mammal. Biophys. J. 11:886-914. 



Cowan, C. C. 1979. Streams: An introduction to fluid flow. An 

instructional module: Physical processes in terrestrial and aquatic 
ecosystems. Ctr. Quantitative Sci., Univ. Washington, Seattle, 

Dawson, W. R. , and R. Templeton. 1966. Physiological responses to 
temperature ir .,h lizard, Crotaphytus eollaris. Physiol. Zool. 
36:219-236. 

Derby, R. W. , and D. M. Gates. 1966. The temperature of tree trunks 
calculated and observed. Amer. J. Bot. 53:580-587. 

Erkine, D. J., and J. R. Spotila. 1977, Heat-energy-budget analysis and 
heat transfer in the largemouth blackbass (Micropterus sal^o^-des) » 
Physiol. Zool. 50:157-16?. 



70 



65 



Gates, D. M. 1968. Transpiration and leaf temperature. Annu. Rev. Plant 
Physiol. 19:211-238. 

Gates, D. M. 1977. Transpiration and leaf temperature. An instructional 
module: Physical processes in terrestrial and aquatic ecos-ystems. 
Ctr. Quantitative Sci., Univ. Washinijton, Seattle. 

■Gates, D. M. and R. D. Stevenson. 1979. Radiation incident on an 

organism. An instructional module: Physical processes in 

terrestrial and aquatic ecosystems. Ctr. Quantitative Sci., Univ. 
Washington, Seattle. 

Gates, D. M. , and C. M. Benedict. 1963. Convection phenomena from plants 
in still air. Amer. J. BoL. 50:563-573. 

Gatesby, R. M. 1977. Conduction nf heat from sheep to ground. Agri. 
Meteor. 18:387-400. ^ s- p ^• 

Golley, F. B. , and J. B. Gentry. 1964. Bioenergetics of the Southern 
harvester ant, Pogonomyrmex badius. Ecology 45:217-225. 

Grace, J. 1977. Plant response to wind. Acad. Press, New York. 204 pp. 

Jackson, W. B. 1957. Microclimatic patterns in the army ant bivouac. 
Ecology 38:276-285. 

Kowalski, G. J. 1978. An analytical and experimental investigation 

of the heat loss through animal fur. Ph.D. dissertation. University 
.of Wisconsin. Madison. 

Kowalski, G. J., and J. W. Mitchell. 1976. Heat transfer from spheres in 
the naturally turbulent, outdoor environment. J. Heat Trans. 
98:649-653. 

Kreith, F. 1973. Principles of heat transfer. Intext Press, New York. 
656 pp. 

Lowry, w. P. 1969. Weather and life. Acad. Press, New York. 305 pp. 

Lueke, R. H., V. Natarajan, and F. E. South. 1^75. A mathematical 
bioth^rmal model of the California sea lion. J. Thermal Biol. 
1:35-34. 

Monteith, J. L. 1973. Pri ciples of environmental physics. Anerican- 
Elsevier Publ. Co., New York. 241 pp. 

Morhardt, S. A., and D. M. Gates. 1974. Energy-exchange analysis of the 
bedding ground squirrel and its habitat. Ecol. Monogr. 44:17-44. 

Mount, L. E. 1968. The climatic physiology of the pig. Edward Arnold, 
London . 

Nobel, P. S. 1975. Effective thickness and resistance of the sir bound- 
ary layer adjacent to spherical plant parts. J. Exp. Bot. 26:120-130. 



66 



Porter, W. P. 1967. Solar radiation through the living body wall of 
vertebrates with emphasis on desert reptiles. Ecol. Monoer. 37: 
273-296. 



^Porter, W. P. 1969. Thermal radiation on metabolic c'.^ambers. T-5cience 
166:115-117. 

•Robinson, N., ed. 1966. Solar radiation. Elsevier Publ. Co., New York. 
347 pp. 

Roseman, L. 1978. Radiation. An instructional module: Physical 

processes in terrestrial and aquatic ecosystems. Ctr. Quantitative 
Sci., Univ. Washington, Seattle. 

Scherba, G. 1962. Mound temperatures of the ant Formioa ulkei (Emery). 
Amer. Mid. Nat. 67:313-385. 

Schraidt-Nielsen , K. 1975. Animal physiology: Adaptation and 
environment* Cambridge Univ. Press, London. 699 pp. 

Schneirla, T. C, R. C. Brown, and F. C. Brown. 1954. The bivouac or 

temporary nest as an adaptive factor in certain terrestrial species 
of army ants. Ecol. Monogr. 24:269-296. 

Siegel, R., and J. R. Howell. 1972. Thermal radiation heat transfer. 
McGraw-Hill Book Co., New York. 

Simpson, J. 1977. Soil heat flow. An instructional module: Physical 
processes in terrestrial and aquatic ecosystems. Ctr. Quantitative 
Sci . , Univ. Washington , Seattle . 

Smart, R. E., and T. R. Sinclair. 1976. Solar heating of grape berries 
and other spherical fruits. Agri . Meteor. 17:241-259. 

Snyder, G. K. , and W. W. Weathers. 1975. Temperature adaptations in 
amphibians. Amer. Nat. 109:93-101. 

Spotlla, J. R., 0. H. Soule, and D. M. Gates. 1972. The biophysical 

ecology of the alligator: Heat energy budgets and climate spaces. 
Ecology 53:1094-110?. 

Stevenson, R. D. 1977a. The First Law of Thermodynamics for ecosystems. 
An instructional module: Physical processes in terrestrial and 
aquatic ecosystems. Ctr. Quantitative Sci., Univ. Washington, 
Seattle. 

Stevenson, R. D. 1977b. Applications of the First Law to ecological 

systems. An instructional module: Physical processes in terrestrial 
and aquatic ecosystems. Ctr. Quantitative Sci., Univ. Washington, 
Seattle. 

Tennekes, H., and J. L. Lumley. 1972. A first course in turbulence. The 
MIT Press, Cambridge, Mass. 300 pp. 



72 



67 



Thorkelson, J., and R. K. Maxwell. 1974. Design and testing of a heat 
transfer model of a raccoon (Procyon lotor) in a closed tree den. 
Ecology 55:29-39. 

Tibbals, E. C. , E. K. Carr, D. M. Gates, and F. Kreith. 1964. "^Radiation 
t and convection in conifers. Amer. J. Bot. 51:529-538. 

-Tracy, C. R. 1976. A model of the dynamic exchanges of water and enerrv 
between a terrestrial amphibian and its environment. Ecol. Monog-. 
46:293-326. 

Van Wijk, W. R. , ed. 1966. Physics of plant environment. North- 
Holland , Amsterdam. 

Welsh, W. R. , and C. R. Tracy. 1977. Respiratory water loss: A 
predictive model. J. Theoret. Biol. 65:253-265. 

Zemansky, M. W. , and H. C. Van Ness. 1966. Basic engineering 
thermodynamics. McGraw-Hill Inc., New York. 380 pp. 




APPENDIX I. Symbols, Units and Dimensions 



'1 



Quantity 



Area in contact with substrate 



Unit 



m 



Dlmenslc 
.2 



Area of emittance 

Percent of total surface area 
exposed to radiation source i 

Constant 

Absorptivity at wavelength A 

Average absorptivity to longwave 
radiation 



in 



Br 
b 
C 

c 



ERIC 



Average absorptivity to shortwave 
radiation 

Breathing rate 

Constant 

Convection 



Heat capacity 
Specific heat 

Water vapor concentration of the 
air 



Water vapor concentration of 
exhaled air 



rain 



-1 



W ra 



-z 



o.-l 



J "C 



kg m 



-1 



kg m 



-1 



74 



-1 



-2 - 

H L ' T 

H M"-^ 0" 



-1 



M L 



-1 



I .1 



►4 



APPENDIX I. Symbols, Units and Dimensions (continued) 

Quanti ty UnU Dimension 

Speed of light 3 x 10^ m s ^ L T ^ 

C Water vapor concentration at the kg ^ ^ ^ ^ 

° surface of the organism 

-3 0 -1 ."^ 
Volumetric specific heat • J m °C H L 0 

2 

c Constant = 2TThc 

1 Li 

c^ Constant - hc^ 

D Characteristic length " m ^ 

DD Degree days °C day 0 T 

E Evaporation M T ^ 

Radiation flux per unit frequency J m » 

E 



, . ' T -3 -1 ht"^t"^ 

Radiation flux per unit wavelength J m s n l i 

Vapor pressure kgm^s^ MLT 



Energy in a photon 



J H 

-1 .-1 

f Frequency s i ^ 

.. G . Conduction W H.T 

-2 -2 

g Acceleration of gravity ni s L T 

Gr Crashof number 



ERIC 
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APPENDIX I. Symbols, Units and Dimensions (continued) 

Quantity Unit 

Coivection coefficient W K 

-34 

Planck's constant 6.63 x 10 J s 

Thermal conductivity W m ^ 

f ""23 ~1 

Boltzmann's constant 1.38 x 10 J K 



Thermal conductance W C 

-2 -1 -k 

Constant for leaf convection 9.1A J ra °C S 

-k -1 

Constant for evaporation flux in 183 S m 
leaf 

Latent heat of evaporation J kg ^ 

Length m 

Length m 

Mass kg 

Metabolism W m 

Molecular weight g mole ^ 
NusseU number 

Wave number m 

Heat energy W m 
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APPENDIX I. Symbols, Units and Dimensions (continued) 



I .1 



Quaiitit 



Absorbed radiation 



Unit 



W 



Dimension 



H T 



-1 



Emitted radiation 



W m 



-2 



-2 -1 

H L T 



Gas constant. 



J mof k"^ 



0 



-1 



Atmospheric radiation (longwave) 

Ground radiation (longwave) 

Shortwave radiation frora the sun 
and sky 

Respiration rate 



W m 



-2 



W m 



-2 



-2 



W m 



3 -1 

m s 



Reynolds number 



Relative humidity 



Resistance 



W"^ m^ °C 



-1 2 

H L T 0 



Reflectivity of the underlying 
surface 



Resistance to water vapor loss 



Internal resistance of leaf 



s m 



-1 



s m 



-1 



T L 



-1 



T L 



-1 



Reflectivity at wavelength A 
Shortwave radiation from the sun 



Fihortwave radiation from the sky 



W m 



-2 



W m 



-2 



-2 -1 
H L T 
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B 



APPENDIX I. Symbols, Units and Dimensions (continued) 
Quantit y Unit 
Temperature . °C 
Boundary temperature "C 



T Air temperature C 

a 



T Maximum daily air temperature C 

amax 



T Minimum daily air temperature C 

amm 



T Body temperature C 

b 



Fruit temperature C 



T Ground temperature C 

g 



Leaf temperature C 



Op 

T Surface temperature v. 
s 

Threshold temperature °C 

T Wall temperature "C 
w 

t Time s 
t(X) Transmittance at wavelength A 

U Internal energy J 

-1 

V • Fluid velocity t^Q m s 

• 3 

Volume m 



APPENDIX I. Symbols, Units and Dimensions (continued) 

Q uantity Unit 

3 

V. Volume of the lungs m 

L 

W Leaf dimension transverse to the m 

wind 

X .V Boundary layer thickness m 

2 Depth ra 

e Emittance 

2 -1 

K Thermal diffusivity in s 

\ Wavelength m 

A Wavelength of maximum radiation ra 

max 

2 -1 

V Kinematic viscosity m s 

p Density Kg m 

-8 -2 -A 

0 Stefan-Boltzmann constant 5.67 x 10 W m K 

, -1 -1 

y Fluid viscosity Kg m s 

M ^ Mass 



0 = Temperature 



